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ABSTRACT 

We have laid down the formahsm and techniques necessary for computing 
the multipole components in a spherical harmonic expansion for bursting 
sources given any specific power spectrum of density perturbations. Using this 
formalism we have explicitly computed and tabulated the expected first few 
multipole components for the GRB distribution using a Cold Dark Matter power 
spectrum. Unfortunately, our analysis leads us to expect that an anisotropy 
signal for this model of structure formation will be below the shot noise level 
for the foreseeable future. Earlier studies by others had claimed that it may be 
possible within a decade or so to probe structure formation using the angular 
distribution of GRBs. We find that while it may be possible to probe the dipole 
due to our motion with respect to the cosmic rest-frame, it does not seem 
feasible to probe the multipole components due to intrinsic fluctuations if the 
actual power spectrum is fixed to within an order of magnitude by the Cold 
Dark Matter model of structure formation. 
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1. Introduction 

The isotropy and paucity of weak bursts suggest that Gamma-ray bursts (GRBs) are 
most hkely cosmological. Analysis of the number vs. peak flux distribution (Cohen & 
Piran, 1995; Fenimore et al., 1993) show that the bursts originate from z 1. This suggests 
that GRBs could be applied as a tool for studying cosmology and the large scale structure 
of the Universe. If GRBs are cosmological, they constitute a homogeneous population with 
a narrow luminosity function (the peak luminosity of GRBs varies by less than factor of 10 
(Cohen & Piran, 1995; Horack & Emshe, 1994; Horack, Emslie & Meegan, 1994; Ulmer, 
Wijers & Fenimore, 1995; Ulmer & Wijers, 1995) that is located at relatively high redshifts 
(Piran, 1992; Mao, & Paczyhski, 1992; Wickramasinghe et al, 1993; Petrosian, 1993; Cohen 
& Piran, 1995). The universe and our Galaxy are transparent to MeV range 7-rays (see e.g. 
Zdziarski & Svensson, 1989). Hence GRBs constitute a unique uniform population which 
does not suffer from any angular distortion due to absorption by the galaxy or by any 
other object. Additionally, GRBs are sampled almost uniformly over 47r steradians. The 
information on matter fluctuations in the Universe is divided to local information obtained 
from red-shift surveys that span distances up to 2; ~ 0.05 or less and to information from 
the CMBR fluctuations that originate at 2; ~ 1000. The recent flnding of KoUat & Piran 
(1996) that GRBs are correlated with Abell clusters suggest that GRBs do follow the matter 
distribution and hence they could provide a unique information about the distribution of 
matter at 2; ~ 1 and on scales that, at present, cannot be explored in any other way. 

Recently, Lahav, Piran & Treyer (1996) developed a scheme for estimating the expected 
spherical harmonics that describe the fluctuations in the intensity of the x-ray background. 
We modify this scheme for estimating the spherical harmonics in the number of bursting 
sources and apply the modifled scheme to the GRB distribution. There are two sources for 
such fluctuations: fluctuations in the number of sources that arises from fluctuations in the 
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matter density or from other effects such as the Sachs- Wolfe effect or the Compton- Getting 
effect due to the motion of the earth and random fluctuations that arises from Poisson 
noise. We estimate the expected fluctuations in the GRBs' angular distribution that arise 
from both sources and we compare these theoretical estimates with recent and future 
observations. The outline of this paper is the following: In section 2. we develop the general 
formalism for estimating the spherical harmonics in a general population of bursting sources 
whose density fluctuates proportionally to the fluctuation in the matter density. Using this 
formalism we estimate in section 3. the number of GRBs needed in order that the real 
angular deviations of the GRB population will be larger than the Poisson noise. In section 
4. we estimate the fluctuations due to the Sachs- Wolfe effect and the Compton-Getting 
effect. In section 5. we compare our results with previous works and with observational 
estimates. 



2. The Formalism 

We consider a general population of bursting sources that originate at cosmological 
distances and follow the matter distribution. The formalism is general but since we apply 
it to gamma-ray bursts (GRBs) we call the sources for simplicity GRBs in the rest of the 
discussion. We expand the number counts of GRBs in a given direction f over the sky, a{r) 
in spherical harmonic expansion: 

(^{r) = J2(^imYim{r) (1) 

Im 

where 

aim^ E (2) 

sources 

To compare the observed distribution with a theoretical one We first estimate the 
number of observed bursts in a given direction A/'(f). To do so we must made several 
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assumption on the nature of the distribution of GRBs. We first assume that there is a 
hnear biasing between the density of GRB sources and the mass fluctuation 

Sx{rc,r) = 6{rc,r), (3) 

where is the bias factor between the density perturbations and the GRB sources, 
6{r) = 6p/p{r) is the density perturbation and Vc is the comoving distance. We also assume 
that the rate of bursts per comoving volume per unit proper time evolves with redshift as 
(1 + zy. At present p is unknown (see e.g. Cohen & Piran, 1995). Finally, we must consider 
the detection procedure of the specific detector. BATSE, for example, is sensitive to C the 
number of photons that arrive on a collecting area A, within a given period of time, AT, in 
a given energy interval centered around E: 

L AAT 

47rr2(l + z)i+" 

where a is the slope of the spectrum. 

In general there exists a sensitivity function, that describes the detectability 

of a burst with a count rate C. Generally, $(C) = 1 for C ^ Cmin and $(C) = for 
C ^ Cmin, where Cmin is the critical sensitivity of the detector. We will usually assume for 
simplicity that the detector behaves like a Heaviside function and $(C) = 1 for C > Cmin 
and $(C) = for C < Cmin- In this case there will be a maximal comoving distance Vc^max 
from which a burst is detected for which C{rc^max) = Cmin- 

Using this assumptions we find that: 

Af{r) = dL^{L) J rldr,n{l+z)'^P~^\l+b-y6{r„r)]<i>[C{L,rc)], (5) 

where Tc is the comoving distance, is the GRB luminosity function, n is an overall 

normalization such that for an observation period Tots, "".(1 + zY /Tobs equals the effective 
density (per unit comoving volume per unit observer time) of bursts. Eq. ^ differs from 
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the corresponding equation for the x-ray background (Lahav, Piran & Treyer, 1996) in two 
ways. First the x-ray background equation, in which we estimate the fluctuations in the 
flux, is weighted by the flux from a given source. There is no such weight here since we 
are interested in numbers of sources. Second, there is additional factor of (1 -|- 2;)"^ in this 
equation since we are interested in bursting sources, while for the x-ray background we were 
considering steady state sources. The predicted harmonic coefficients are: 

aim = j drM{v)YCm{i) = 
j dvj^ dL^L) J rldrc n{l + z)^p-'^[1 + b,S{r„ f)]<l> [re,^„,(L) - r,] F4(f), (6) 

If N is the total number of observed bursts, we can use the monopole (/ = 0), with 
^00 = (47r)~^/^, to fix the normalization factor n in the following way: 

TV ^ (47r)V2^oo- (7) 

For the special case of standard candles, $(!/) = S{L — Lq), with no source evolution (p = 0) 
we have 



N = 



3 

c.max 



1 + 



(8) 



20c2 4c 

where Hq is the Hubble constant and Tcmax is the maximal comoving distance from which a 
GRB with luminosity Lq can be detected. Thus, for an Einstein-de Sitter universe we have: 

3N /Hn\^ 1 1 



n 



327r V c 7 \i 1 



a 



0.1 + 7T3^ + 



0.6 I 0.3 



(9) 



where Zmax{L) = z{rc,max{L)) is the maximal redshift from which a GRB with luminosity L 
can be observed. 

The fluctuations in the background are expressed by higher harmonics I > 0: 

aim^nj dL^L) j dv j r^dr, $[C(L, rJlMl^c ?) Y;:^{i){l + zf-\ (10) 
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Expressing the fluctuations in the matter density in terms of the Fourier components Sk 
(see Eq. H) and using the orthogonality of Spherical Harmonics / dujYim{r)YiT^^,(r) = b"^!^' 
we find: 

aim = ^(^')*^ / rf'A;6^5k(^)l^™(k) J dL^L) J r^dr, <!>[C{L,r,)]ji{kr,){l + zy-\ (11) 

We are interested in sources at 2; ~ 1 and at large wavelengths. The relevant 
fiuctuations that contribute to this integral are still in the linear regime. This means that: 

Si,{z)=6^o/{l + z) (12) 

where ^ko is the present linear amplitude of fluctuations normalized to fit the CMB 



observations. Eq. |T2| is gauge dependent (in some gauges perturbation that are larger that 
the horizon do not grow). However, it can be checked numerically that the contribution of 
very large scale perturbations is small and hence this gauge choice is not important. 

Taking the mean-square values and using Parseval's theorem and 

{5^{z)5Uz)) =(27r)3p(A;,^)5(3)(k-k') (13) 

we obtain 

{\aim?) = - b',i—fn' / dkePik)\<ili{k)\\ (14) 
The function \l'i(/c) is defined (for an Einstein -de Sitter universe) as: 



^i{k) = I dm{L) j dz{i + zf-^i^ 



1 ' 



jKA;r,)$[C(L,z)]. (15) 



3. Estimates for GRBs 

We consider standard candles GRBs with no density evolution i.e. p = 0. For simplicity 
we use in this section $(C) = 1 for (7 > Cmin and $(C) = otherwise. The power-spectrum 
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for CDM can be written down in the form (Padmanabhan, 1993): 

^^^^ ^ {l + Bk + Ck^/^ + Dk^y ^^^^ 

For a scale-invariant or Harrison-Zel'dovich spectrum, n = 1. Further, by fitting to available 
data we get the best-fit values of the parameters as follows: A = (24Mpc)^, B = 1.7Mpc, 
C = dMpc^^"^ and D = IMpc^. Let us scale out the dimensions of k in the following way, 

k = k{100Mpc) (17) 

By doing this rescaling the expression for (la^mP) can be written completely in terms of 
dimensionless numbers in the following form: 

(3.0 X 10-5)A^2 ,2 /■ .r^^l^K'^)!^ no^ 



= — ^ / dk 



Zmax 



where 

D(k) = (1 + 0.017fc + 0.009p/2 + O.OOOlP)^. (19) 
Thus, the quantity which needs to be estimated numerically is: 

/(Z, Zma.) ^ I dk^^^^Ml (20) 

Let us now obtain an order of magnitude estimate of the above effect and see whether 
this is going to be detectable in the near future (Tegmark et al, 1995; 1996). Let us re-write 



the earlier equation |18| in a form convenient for tabulation of values for various values of 



(lOimP) = PiZmax)N'^ bp{l, Zmax) (21) 

Here p{zmax) is the numerical pre-factor which depends on Zmax, which we shall tabulate 
below. 
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TABLE 1 



^max 


0.2 


0.5 


1.0 


2.0 


P^^max) 


9.0 X 10^ 


1.4 


1.3 X 10-1 


2.4 X 10-2 



One can estimate the magnitude of /(/, z^ax) for various values of I. 



TABLE 2 





0.2 


0.5 


1.0 


2.0 




3.09 X 10-9 


8.20 X 10-9 


1.064 X 10-8 


1.036 X 10-8 


-'- V — - -iiKi.r j 


3.90 X IQ-" 


1.061 X IQ-^ 


1.139 X IQ-^ 


1.187 X 10-8 


-^(3, Zmax') 


4.72 X lO-'' 


1.310 X 10-^ 


1.826 X 10-8 


1.949 X 10-8 




5.51 X 10-9 


1.559 X 10-^ 


2.214 X 10-8 


2.411 X 10-8 


-^(5, Zjyidx) 


6.29 X 10-9 


1.804 X 10-^ 


2.598 X 10-8 


2.868 X 10-8 



Wc, of course have to compare the magnitude of the signal given in the previous 
equation to the shot noise. The magnitude of the shot noise in our case is given by : 

\ r fi TV 

<|a,„R.„ = -/<iK^=j^ (22) 

where N is the number of bursts observed. Let us now estimate what is the number of 
bursts one needs to get to observe a signal above noise. This means we want, 

p{Zmax)N^ hp{l, Zmax) > ^ (23) 

We can obtain an estimate for the number of bursts required to get a signal above 
noise for any desired assumption about the location of the most distant sources of GRBs. 
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in particular let us obtain the required number of bursts if we assume the most distant 
sources of GRBs is at z^ax = 0.2. We start at Eq. ^ and use the value of p{zmax) given in 
Table 1 to get the relation, 

ilaiml'') = iV'(9.0 X 10^) 0.2) (24) 

One can numerically compute the magnitude of /(/, 0.2) for various values of I. We did 
this for 1 = 1 to 1 = 5. The magnitude does not change significantly and is given by, 
/ ~ 4 X 10^^. Thus the relevant inequality now becomes, 

N\9.0xl0^)bp{l,0.2)> ^ (25) 

Thus if the most distant GRBs are at Zmax = 0.2 then the number of bursts required to get 
a signal over noise must be given by > 2 x 10^. In general, if the most distant GRBs 
are at Zmax then the number of bursts required to get a signal over noise must be given by 
> Nc{l, Zmax)b~'^, where 



(26) 



Below we tabulate the Number of bursts, Nc{l, Zmax) required to observe the multipole 
component / for a given maximum redshift of sources Zmax'- 
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'^max 


0.2 


0.5 


1.0 


2.0 


A r / -1 \ 


2.9 X 10^ 


6.9 X 10^ 


5.8 X 10'' 


3.2 X 10^ 


■^ci^i ^max^ 


2.3 X 10^ 


5.4 X 10^ 


4.3 X 10^ 


2.2 X 10^ 


-^c(J^j ^max) 


1.9 X 10^ 


4.3 X 10^ 


3.4 X 10^ 


1.7 X 10^ 


-^c(45 ZjYiax) 


1.6 X 10^ 


3.6 X 10^ 


2.8 X 10^ 


1.4 X 10^ 


-^c(5, ZfYiax^ 


1.4 X 10^ 


3.2 X 10*^ 


2.4 X 10^ 


1.2 X 10^ 



These results are perhaps discouraging because it means that the expected signal is not 
detectable in the foreseeable future. However, given the fact that there are ongoing efforts 
to measure the multipole components of the GRB distributions (Tegmark et al, 1995;1996), 
this is still a significant result. It means that if a signal is detected it is not due to the effect 
discussed here by us and an alternative explanation must then be sought for the signal. 

4. Sachs- Wolfe and Compton-Getting 

We turn now to fluctuations in the redshift that influence the observed count rate and 
induce fluctuations in the number density given by Eq. (|5|). So far we have introduced 
n = no + Sn where uq is the background density and 6n are the fluctuations and we were 
interested in those fluctuations (we assumed that 6n = b^yUoSp where is a biasing factor 
and 6p the fluctuations in the matter density. Now we consider fluctuations in the count 
rate C and in the redshift factor {1 + z) that arise due to Doppler motions of the sources 
and the observer (Vq — f||) (where Vq is our motion and v\\ is the motion of the source), 
or due to fluctuations in the gravitational potential S(j). For simplicity we discuss in this 
section a standard candle population. The generalization to a luminosity function is trivial. 

The variation in the factor of {1 + z) causes fluctuations in the time dilation between 
the sources and the observer and hence fluctuations in the observed rate of events from 
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different regions of space-time. These lead to the following fluctuations in the angular count 
rate: 



The first term Vq is a constant vector that gives rise, just like in the CMBR, to the Dipole 

and a second order quadrupolc term: 



Note that for sources that are at 2; ~ 1 the intrinsic fluctuations discussed earlier (Eq. (21) 
and Tables 1 and 2) give a dipole at the 10"^ level. The Doppler dipole discussed above 
in Eq. (28) is at a few times 10"^ level. Thus the intrinsic dipole is comparable to but 
about an order of magnitude below the Doppler dipole. However, by examining the same 
equations and tables as above we see that the intrinsic quadrupole is clearly much greater 
than the Doppler quadrupole. In general the higher Doppler multipole components are 
suppressed by increasing powers of (Kfts/c) whereas the higher multipole moments of the 
intrinsic fluctuations are of the same order as the lower ones (see Table 2). 

The second and third terms, f|| and Sip, reflect the fluctuations of the velocity fleld and 
the gravitational field. Like in the discussion of the density fiuctuations we turn to Fourier 
space. Using: 




(27) 




(28) 




(29) 



for Sep and 





k 



(30) 



for the line-of-sight velocity we obtain: 




(31) 



and (Fisher, Scharf & Lahav, 1994) 




(32) 
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Substitution of Eqs. ^ and |3^ into ^ yields: 



(33) 



Im 



To calculate the corresponding aim we expand 5A/'(f), expressed in Eq. ^ : 



O'lmA 



2n^ 



n d k F;„(k)5ko 



^0 VI/ M ^°vl/ 



where 



and 



-no "'0 



p-5/2 



n 2 



Hq Jo 



1 



ii{krc) , 



fiikrc) ■ 



(34) 



(35) 



(36) 



The fluctuations in the gravitational fleld and in the velocity fleld also influence the 
detectability of bursts that are marginally detectable and this leads to another effect. The 
BATSE detector is sensitive to C the number of photons in a given energy interval that 
arrive to the detector in a given period of time. The count rate C varies due to variation in 
z which arise from the Doppler effect and the gravitational redshift. Using Eq. ^ we flnd: 

dC (1 + a)C 



dz 



(37) 



Substitution of Eq. 37 to Eq. ^ and linearizing we obtain: 



SAfo 



nil + z 



dC 



+ a)C[-(K) - t'||)/c - m-icyidr, 



(38) 



For simplicity we assume that our detector has a sharp cutoff at Cmin- Generalization 
to a smooth cutoff it trivial. In this case $ is a Heaviside function, |^ = 6{C — Cmin) and 
the integral becomes a surface term! 
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We can express the volume element as a Jacobian times dC. Thus, we arrive at the 
following trivial integral for 6^/2: 



p \ roo IT) p 
5M2={ — ] / ——MC-Crmn)\{Vo-V\\)/c + 5ct>/^c' 

Mq/ Jo [I + z) l 



1 



dC 



n 



{Vo ))/c + 6(f){Zmax)/^C 



Zmax) and r^{zmax) are all calculated on the surface Zmax which 



Mo J {l + Zmax^/^ 

The quantities: v\\{zmax), 
corresponds to Cmin with our standard candle luminosity. 

The contribution from Vq is again a dipole and a second order quadrupole term: 
ANd2 (Vobs/c) cos{e) + {Vobs/cy[l - cos2(^)/2] 



(39) 



N 



a/1 + 



0.1(1 + Z,nax) + 0.3a/1 + Zmax + 0.6 



(40) 



We discuss the magnitude of these terms as well as those given in Eq. ^ in the next section. 

The fluctuation terms can be dealt in the same manner that we have used earlier. 
First, we substitute Eqs. ^ and |32| into Eq. Then, after multiplying by Yim{r) and 
integrating over dQ we get: 



27r2 



n / d-^k Y(^(k)(5ko 



Zk'^c^ kc 



where: 



and 



^M-,5(A:)=4(— )3(l + ^„,,)f-3/2 
-no 



1 - 



= 4(— )=^(1 + z 



Zma 
1 ^2 



ii{krc 



j[{krc 



(41) 



(42) 



(43) 



We can combine now all fluctuating terms (Eqs. 0, |3^ and ^) to a single equation for 
aim which resembles Eq. |ll] with ^f; replaced by + ^!i^sw + '^i,sw,s + ^i,v + '^i,v,s- 



air, 



27r2 



n / d^k Yim(k)5ko 



H"^ H 



2k^C^ 



(44) 
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Now, we continue and calculate {aimO*im) taking the mean-square values and using Parseval's 
theorem and Eq. O: 



J^2 jj 



(45) 



^From the above expressions we can calculate the most general r.m.s. amplitude 
(aimCi'im)- It is worthwhile, however, to estimate the order of magnitude of the relevant 
terms. To do so we recall that if A is the length scale being probed then (Padmanabhan, 
1993): 

and 

(..-„,„/e.g)^(iLA). 

These estimates suggest that the Sachs- Wolfe terms in Eq. ^ are smaller than the intrinsic 
fluctuation term by a factor of {H\/cY. The velocity terms are smaller than the intrinsic 
fluctuation by a factor of H\/c. 

Comparison of the surface terms in the Sachs- Wolfe and the velocity terms to the 
volume integral show that these terms are comparable for small Zmax and become smaller 
as Zmax becomes larger. This is intuitively explained by the fact that at lower Zmax there 
are relatively more faint bursts than at higher z^ax- 

^From our expressions it is clear that there is a systematic expansion in terms of the 
parameter [Ho/kc). These results, show that the Sachs- Wolfe and the velocity terms can be 
ignored if one is interested in the lowest order effect and an order of magnitude calculation 
as we have been interested here. The difference between GRBs and CMBR is that the first 
arise from sources at relatively low red- shift, for which the intrinsic fluctuations have grown 
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and are much more important than the fluctuations induced by a varying gravitational fleld 
or from the random pecuhar velocity of the sources. 

5. Discussion 

We have laid down the formalism and techniques for evaluating the expected multipole 
components of the GRB count distribution. The largest term, which is a dipole of 
order 10~^, almost two orders of magnitude larger then all other terms arises due to the 
Compton- Getting effect resulting from our motion relative to the GRB distribution. This 
term is not new. It was first calculated by Maoz (1994) who estimated to be of order 
Ri 10"^. Later Scharf et al. (1995) repeated the calculations and have considered also the 
fluence weighted dipole (which we don't consider here). 

Smaller multipole components arise from intrinsic fluctuations of the source 
distribution that are related to the power spectrum of density perturbations. We have 
explicitly computed the first few multipole components for a Cold Dark Matter power 
spectrum. These are in fact the very multipole components that are least likely to be 
affected by the positional inaccuracy in locating GRBs. We feel that the expression of the 
fluctuations in the GRB distribution in terms of a multipole moments is particularly useful 
in view of the general usage of these moments to express fluctuations in the CMBR (Bond 
& Efstathiou, 1987) and other background flelds (Peebles, 1973; Lahav, 1994; Fisher, 1993; 
Heavens & Taylor, 1994; Lahav, Piran & Treyer, 1996). There have now been serious and 
exhaustive efforts (Tegmark et al, 1995; 1996) to examine the deviation from isotropy of 
GRBs using an expansion in terms of spherical harmonics. Further, this observational 
program is likely to continue into the future and so a calculation of the expected signal in 
terms of various multipole components is a useful calculation to do. 
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Our work would certainly be incomplete if we failed to mention the earlier work of 
Lamb & Quashnock (1993) who were the first to address the issue of whether it is possible 
to probe Large Scale Structure using the distribution of Gamma Ray Bursts, and did 
not compare our results to their analysis and conclusions. LQ use the angular two-point 
correlation function to examine the expected deviations from isotropy, whereas we have 
used an expansion in terms of spherical harmonics. While undoubtedly both formalisms 
have their usefulness we have already mentioned that we do feel that an expansion in terms 
of spherical harmonics is better suited to address the issue of deviations from isotropy. 
However, having said all that in favor of doing a spherical harmonic expansion of the GRB 
distribution it still remains a true statement that such conclusions as the number of bursts 
required to see an anisotropy signal over the noise should be basis independent. 

We now turn to comparing our results about the number of bursts required to see 
a signal above noise to those obtained by LQ. The central observational quantity in the 
formalism that LQ use is the mean of the angular two-point correlation function, w{9). 
The statistical uncertainty in w{6) is denoted by 6w{6). In order to detect a signal it is 
therefore necessary to require Sw{9)/w{9) < 1. The above inequality can be transformed 
into an inequality which determines the minimum number of bursts required to detect a 
signal due to the actual anisotropy of the distribution as opposed to the expected statistical 
fluctuations. This bound is written down by LQ as 



where Nf, is the number of bursts, D is the distance to the edge of the source distribution, 
ko is the largest wavenumber that can be probed, P{k) is the power spectrum of density 
perturbations, C is a constant introduced by LQ and I{ko) is the following integral 
introduced by LQ, 



Nb> 




(48) 



CI{ko)P{ko) ' 




(49) 
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Here W{k) is a window function introduced by LQ. While LQ are somewhat ambiguous 
about the exact numerical value of the window function used to arrive at their results, they 
do argue that in the regime < k < [D9)~^ the window function can be computed 
analytically tending for large values of kD to Further, they argue that most of the 
contribution to the integral, I{ko) comes from wavenumbers k ~ ko- Using these inputs we 
obtain an order of magnitude estimate for I{ko) to be 



We also need to input the power spectrum to extract the bound on the number of 
bursts for a detectable signal. To make a fair comparison let us do this along the hues 
suggested by Figure 1 of LQ. Thus we take, ko = 0.1hMpc~^ and P{ko) ~ 5 x lO^h^Mpc^. 
LQ also argue that C ~ 1.5 in order for their convolution to match those from Limber 
equation for small values of 9. Further they choose D — lh~^Gpc for illustrative purposes. 
Note that this corresponds to Zmax — 0.3. Putting in all these values we get the bound on 
the number of bursts, for a detectable signal to be. 



This result agrees in order of magnitude with our results for comparable Zmax- Thus we 
see that within the context of our analysis both the results of the LQ formalism and the 
spherical harmonic expansion agree. However, there is one claim made by LQ with which 
we disagree. Namely, they claim that if an anisotropy signal is not detected with ~ 3000 
bursts then it would cast a doubt on the cosmological origin of GRBs. Such a statement we 
think is unwarranted. Certainly if one stays within the context of the Cold Dark Matter 
models whose spectrum we have used in our analysis and the maximum redshift from which 
GRBs originate Zmax > 0.2 there would be no reason to expect an anisotropy signal above 
that of shot noise with only ~ 3000 bursts. 



I{ko) ~ 5MkoD . 



(50) 



iVh > 2 X 10^ . 



(51) 
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Let us now summarize our results in the light of the observational situation. We first 
note that BATSE has accumulated ~ 1000 bursts in 3 years of operation. If it continues 
operating for a period of 10 years it will accumulate ~ 3000 bursts. In light of our results 
displayed in Table 3 this means that if GRBs originate at redshifts Zmax > 0.2, and the 
Cold Dark Matter power spectrum is right to within an order of magnitude, then we are 
unlikely to see an anisotropy signal in GRB distributions above the shot noise within the 
next 10 years or so. Given that there are now serious and exhaustive efforts to measure and 
analyze the spherical harmonic components of the GRB distribution (Tegmark et al, 1995; 
1996), we think an accurate estimate of the signal expected from Large Scale Structure was 
necessary. 

We have laid down in this paper the formalism and techniques necessary for computing 
the various multipole components in a spherical harmonic expansion for bursting sources 
given any specific power spectrum of density perturbations. We have further explicitly 
computed and tabulated the first few multipole components for GRBs distribution using a 
Cold Dark Matter power spectrum. Unfortunately, our analysis leads us to expect that an 
anisotropy signal for this model of structure formation will be below the shot noise level for 
the foreseeable future. 
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Appendix: Fourier Preliminaries 

Fourier decomposition of the density perturbations 

and tlie inverse: 

(5k = y d^r 5{r) e"^ ' . (2) 

Rayleigh expansion: 

e^^- = 4nJ2^'Jl{kr)Y,U^)Y:M . (3) 

Im 

Combining Eqs. (|l|) and (|[) : 

'^(r) = ^ Ei^TYLi^) I d'k M^(k) Ukr) . (4) 
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